Electron-phonon effects on spin-orbit split bands of two dimensional systems 
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The electronic self-energy is studied for a two dimensional electron gas coupled to a spin-orbit 
Rashba field and interacting with dispersionless phonons. For the case of a momentum independent 
electron-phonon coupling (Holstein model) we solve numerically the self-consistent non-crossing 
approximation for the self-energy and calculate the electron mass enhancement m*/m and the 
spectral properties. We find that, even for nominal weak electron-phonon interaction, for strong 
spin-orbit couplings the electrons behave as effectively strongly coupled to the phonons. We interpret 
this result by a topological change of the Fermi surface occurring at sufficiently strong spin-orbit 
coupling, which induces a square-root divergence in the electronic density of states at low energies. 
We provide results for m* /m and for the density of states of the interacting electrons for several 
values of the electron filling and of the spin-orbit interaction. 

PACS numbers: 71.38.-k, 71.70.Ej, 73.20.At 



I. INTRODUCTION 

Prompted by considerable technological interests, the 
physics of itinerant electrons coupled to spin-orbit (SO) 
potentials has been the subject of extensive investiga- 
tions in recent years In materials of interest, the main 
sources of SO coupling are the Rashba interaction arising 
from structural inversion asymmetries of low-dimensional 
structures^ 3 - and the Dresselhaus interaction present in 
bulk crystals lacking inversion symmetry^ Depending on 
the material characteristics, one of the above interac- 
tions, or even both, may be present, lifting the spin de- 
generacy of the electron dispersion. When measured at 
the Fermi level, the resulting energy splitting, A so , is 
commonly used to estimate the strength of the SO inter- 
action. 

In narrow-gap III-V semiconductor-based heterostruc- 
tures, such as GaAs and InAs quantum wells, A so 
is a few meV, while in II- VI quantum wells A so is 
greatly enhanced. For example the heavy-hole conduc- 
tion band of HgTc displays SO splitting values ranging 
between 10 — 17 meV and 30 meV^ Much stronger 
SO splittings have been observed in the surface states 
of metals^ and semimetals^^, and the corresponding A so 
may be so large, e.g. A so ~ 110 meV in Au(lll))£ 
that the possibility of detecting SO split image states 
has been recently put forward.— Other systems display- 
ing giant SO splittings arc surface alloys as, for exam- 
ple, Li/W(110),iiPb/Ag(lll)^Hand Bi/Ag(lll)^ior 
even onc-dimcnsional structures such as Au chains in vici- 
nal Si(lll) surfaces^ For such low-dimensional or struc- 
tured materials, the SO interaction is of Rashba type, but 
large SO splittings have been found (or predicted) also 
in bulk crystals, where the Dresselhaus interaction leads 
to A so as large as 200 meV in non-centrosymmetric su- 
perconductors CePt 3 Si,i&±l Li 2 Pd 3 B, and Li 2 P7 3 B3^ 



Such strong SO couplings may possibly have interest- 
ing applications in spintronic devices, but represent also a 
compelling and challenging problem from the theoretical 
standpoint, in particular when A so is no longer the small- 
est energy scale in the system, as in III-V semiconductor 
heterostuctures where A so rj 1-5 meV, but competes 
in magnitude with other characteristic energy scales such 
as the phonon frequency or the Fermi energy. From this 
perspective, systems like the Bi/Ag(lll) surface alloy, 
which shows bands split by about 200 meV ) 13 ' 14 are par- 
ticularly promising, given also the alleged possibility of 
tuning, by Pb doping, the Fermi energy Ep to values 
lower than the SO energy splitting^ 

A few novel and interesting features arising from strong 
SO splittings have already been investigated theoretically 
in the literature. For example, in Ref. [2~fl| it has been 
demonstrated that the Rashba SO coupling induces an 
infinite number of bound states in two dimensions, even 
for short ranged impurity potentials, while in a recent 
work we have shown that the superconducting critical 
temperature of a low-density two-dimensional (2D) elec- 
tron gas can be significantly enhanced by the Rashba 
interaction^ 2 - Both phenomena discussed in Refs. [2lH22| 
can be understood in terms of a SO induced topological 
change of the Fermi surface, which gives rise to an effec- 
tive reduction of dimensionality of the electronic density 
of states for Ep sufficiently smaller than the SO charac- 
teristic energy 

In this paper we analyze the effects of such topologi- 
cal change of the Fermi surface on the electron-phonon 
(cl-ph) problem of 2D systems. In particular, we study 
one-particle spectral properties and extract the combined 
el-ph and SO effects on the electronic effective mass m* 
and on the interacting density of states (DOS). We show 
that, even for weak or moderate couplings to phonons, 
the effective reduction of the bare DOS induced by the 
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FIG. 1: (color online), (a): electron dispersion = ^(fci 
fco) 2 for a, spin-orbit split electron gas. The energy Eq = -^i^ 2 
is a measure of the spin-orbit interaction and is equivalent 
to the minimum interband excitation energy for an electron 
sitting at the bottom of the lower band. The upper and lower 
horizontal dashed lines indicate the position of the Fermi level 
for Ef > Eo and Ef < Eo, respectively. Also shown are the 
corresponding Fermi circles with occupied states drawn by 
gray (colored) regions, (b): density of states plotted from 
Eqs.flSl7]). 



Rashba interaction leads to a strong increase of m* , and 
to phonon satellite peaks in the interacting DOS, which 
are typical signatures of an effectively strong el-ph cou- 
pling. Due to the two-dimensionality of our model, and 
to the Rashba type of SO coupling, our results could 
be relevant for both metal and semimetal surface states, 
for which the el-ph interaction has been shown to be 
relevant ) 9 ' 23 i 24 ' 25 i 26 and for surface superconductors ^1 
with the hypothesis that pairing is provided by the cou- 
pling to phonons. 



of the confining potential. Here we consider a linear 
Rashba model for the SO interaction 
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where 7 is the SO coupling constant. Furthermore, we 
assume that the unperturbed electron band is parabolic, 
e k = h 2 k 2 /2m, where m is the band mass of the electron. 
Apart from a constant shift Eq (defined below) which can 
be absorbed in the chemical potential, the eigenvalues of 
Eqs.((T|) and © are: 
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where k = |k|, s = ± is the chirality number, and fco is 
the Rashba momentum 
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The two electron branches are plotted in Fig. HJa) 
in units of the Rashba energy 
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which corresponds to the energy difference between the 
degeneracy point at k = and the bottom of the lower 
band at k = fco. In Fig. [3a) we indicate also the Fermi 
levels for the Ep > Eq and Ep < Eq cases (horizon- 
tal dashed lines) which represent two qualitatively dif- 
ferent situations. For Ep > Eq, the Fermi level crosses 
bands of different chirality and the corresponding Fermi 
sea is given by the area of two concentric Fermi circles, 
as sketched in Fig. |Tfa) . In this case, the corresponding 
DOS for each sub-band is 



N±(E F ) 




for E F > E Q (6) 



II. RASHBA-HOLSTEIN MODEL 

Two-dimensional quantum wells, with strong and 
asymmetric confining potentials, and surface states with 
weak or negligible coupling to the bulk can be satisfacto- 
rily represented by the following 2D electron hamiltonian 
with SO interaction 
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where c ks (c ks ) is the creation (annihilation) operator 
for an electron with momentum k = (k x ,k y ) and spin 
index a =f, J.. In the above expression, e k is the elec- 
tron dispersion in the absence of SO coupling, er is the 
spin-vector operator with components (a x , a y , a z ) given 
by the Pauli matrices, and f2 k is a k dependent SO pseu- 
dopotcntial arising from the asymmetry in the z-dircction 



where No = m/2irh 2 is the DOS per spin direction with 
zero SO coupling. The sum over the two chiral states, 
N(Ep) = N + (E F ) + N-(Ep), is therefore identical to 
the total DOS, 2A*o, of a 2D electron gas without SO 
interaction [Fig. |TJt>)] . Furthermore, in the Ep ^> Eq 
regime, one has N±(Ep) ~ No, and the dispersions of the 
low excitations in the vicinity of Ep can be approximated 
by vp{k — kp)±A so /2, where vp and kp are respectively 
the Fermi velocity and momentum in the absence of SO 
interaction and A so = 2"fkp is the SO energy splitting. 
This is the quantity which is usually used to quantify the 
SO strength in semiconductors such as GaAs and InAs. 

For Ep < Eq the situation is drastically different. In 
this case in fact, as shown in Fig. [IJa), the Fermi level 
crosses only the s = — 1 band but, since E^ has a min- 
imum at k = fco 0, the Fermi surface is still consti- 
tuted by two concentric circles. The resulting Fermi sea 
is therefore given by the area of the annulus comprised by 
the two circles and in the limit of Ep — > 0, with Eq ^ 0, 
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the Fermi surface Sf coalesces into a circle of radius fco, 
Sp = 2irk , while the Fermi velocity vp vanishes as \[Ep~. 
Since N(E F ) on Sf/v f , the resulting DOS is therefore^ 



N{Ep) = N-(Ep) = 2Nq\ -=— for E F <E . (7) 

V E F 

As we shall see in the following, the one-dimensional-likc 
singularity of Eq.© has important and peculiar effects 
on the low-energy properties of the system, in contrast 
with the Ef > Eg case, for which the corresponding DOS 
is featureless. 

Let us introduce now the coupling to the phononic de- 
grees of freedom. In the present paper, we consider the 
following Holstein-type of interaction hamiltonian 
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a Ck'a(4_k' + a k'-k), (8) 



where a q (a q ) is the creation (annihilation) operator for a 
phonon with momentum q, u>q is a dispersionless phonon 
frequency, and g is the momentum independent el-ph ma- 
trix element. As will become clear in the following, the 
choice of the momentum independent quantities luq and 
g is convenient for the calculation of the self-energy, and 
permits a more direct evaluation of the effects of the SO 
interaction on the cl-ph properties. The present analysis 
is therefore a starting point for more general formulations 
of the cl-ph hamiltonian. 

The thermal Green's function of the electrons sub- 
jected to the total hamiltonian H = H a + satisfies 
the following Dyson equation 



G(k,iw„) = [G 1 (k,io; n ) - S(k,iw„)] 



(9) 



where cu n = (2n + 1)ttT is a Fermionic Matsubara fre- 
quency and T is the temperature. Go(k, iuj n ) is the non- 
interacting electron propagator and S(k, iu n ) is the self- 
energy due to the coupling with phonons. Due to the SO 
interaction appearing in H , these quantities are 2x2 
matrices in the spin sub-space. From Eqs. (Q]) and 
the non-interacting propagator is 
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G (k, ioj n ) = ^ ^ (1 + sO k • a)G (E s k ,iu n ), (10) 

1 s=±l 



where Ct^ cr = k x a y — k y a x , and Go(E kl iuj n ) = l/(iuj n — 
E k + fi), where /j, is the chemical potential. 

For the evaluation of the self-energy, we shall consider 
a self-consistent Born approximation (non-crossing ap- 
proximation) which neglects all el-ph vertex corrections. 
Furthermore, we shall not consider many-body correc- 
tions to the phonon propagator. These limitations will 
be discussed in SecfVland. for the moment, it suffices to 
keep in mind that this approximation scheme should be 
not too poor as long as the coupling to the phonons is 
sufficiently weak. Hence, given the phonon propagator 
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the resulting electron self-energy matrix in the non- 
crossing approximation reduces to: 



S(k, ioj n ) 
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D(iuj n - i(j m )G(k',iu) m ), 



(12) 

where A = 2g 2 A /u>o is the el-ph coupling constant. 
From Eq. (|12[) it is clear that, due to the momentum in- 
dependence of the el-ph interaction, the self-energy (fl"2|) 
depends only upon the frequency. Furthermore, by sub- 
stituting Go(k',iw m ) for G(k',iw„ i ) in Eq. (fT2")) . the re- 
sulting second-order self-energy is diagonal in the spin 
space. This holds true for all orders of iteration, so that 
S(k, iuj n ) = ~S(iu) n )l, where 1 is the unit matrix. The 
Green's function (TQl) can therefore be rewritten as 



G(k, ioj n ) 



-y 



where 



G(E s k ,iu n ) 



(l + sf2 k -<r)G(£|,j<), (13) 



(14) 
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is the electron propagator in the chiral basis for the in- 
teracting case while the self-energy is 
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In the above expression, we have introduced an upper 
momentum cut-off k c which prevents the integral over k 
from diverging. Such divergence is an artifact due to the 
use of a momentum independent el-ph matrix element g 
in Eq.© and of the electron gas model of Hq. On physi- 
cal grounds, the introduction of k c is equivalent therefore 
to defining a finite Brillouin zone of area nk 2 or, cquiv- 
alcntly, a finite bandwidth E c = h 2 k 2 /2m when Eq = 0. 
In the following, E c will be chosen to be much larger than 
the other relevant energy scales of the system (E c ^> cj , 
Eo, Ep). A finite k c , or E c , permits also to define a fi- 
nite electron density p e = J2 a J dk/ (27r) 2 (c ] { co .Ck CT ) which, 
relative to the cut-off k c , becomes 
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where o + is an infinitesimal positive quantity and 
the second equality has been obtained by using 
TE n G(EliLo n )e^° + = l/2 + Tj: n ReG(E s k ,iu> n )& 
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FIG. 2: Electron density number function of the bare 

chemical potential /j, for A = 0.5 and several values of the 
SO interaction Eq. The temperature is T = 0.02^o and the 
energy cut-off is E c = lOO^o. In the inset n e is plotted for 
A = and for zero temperature. 

In the following, we shall present results in terms of the 
electron number density 



(18) 

2 (n e = 0) for 



which attains the limiting value n e 
completely filled (empty) bands. 

Before turning to the next sections, where we present 
our numerical results, it is worthwhile showing how the 
SO effects on the DOS enter the self-energy function. By 
transforming the integration over k in an integration over 
the energy, Eg. ([16]) can be rewritten as follows: 



g(iu m ) 



dE Po (E)G(E,iu; m ), 
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where, for simplicity, terms of order \J Eq / ' E c have been 
omitted in the upper limit of integration, and 



Po(E) = 



N S (E) 

2Nn 



1 for E > E, 



for E < E 



(20) 



is the reduced non-interacting DOS obtained from 
Eqs.© and ([7]). From the above expressions, it is there- 
fore straightforward to realize the importance on the el- 
ph properties of the square root singularity of the DOS 
at low energies. As we shall see in the following, the ef- 
fective electron mass and the electron spectral properties 
in the presence of SO interaction will differ qualitatively 
from the corresponding results for Eq = 0. 



III. EFFECTIVE MASS 

The integration over the momenta appearing in 
Eq. (|16[) . or, equivalcntly, the integration over the energy 



FIG. 3: Electron effective mass m* as a function of the elec- 
tron density number n e for A = 0.5 and for several values of 
the SO interaction Eq. 



in Eq.(|19|). can be carried out analytically, leaving only 
the summation over the Matsubara frequency to be per- 
formed numerically. Hence, for fixed values of A, luq and 
Eq, the electron self-energy Y,(iu> n ) is obtained by iter- 
ation of Eqs.([Ti]). (fi~5|) . and (T6]), while Eq.|l7|) is used 
to extract the corresponding electron density for a given 
value of /x. For all cases we have set E c = lOOwo and 
T = 0.02wo, which is low enough to be representative of 
the zero temperature case. In Fig. [2]we report the calcu- 
lated values of n e , Eq. (fT5)) . for A = 0.5 and for different 
values of the SO energy Eq- For comparison, we report 
in the inset of Fig. [5] the corresponding density values 
for A = and at zero temperature. For E = 0, n e de- 
creases almost linearly as \x is reduced, as expected for 
a constant DOS in 2D (see inset), but the zero density 
limit n e = (extracted in the T — ► limit) is reached 
only for fi = /iq ~ — 1.023a->o, which is lower than the 
non-interacting zero-density value jj, = 0. This energy 
decrease represents the ground state energy of a single 
electron in interaction with phonons and provides a mea- 
sure of the strength of the el-ph interaction. For non-zero 
SO coupling, E > 0, two features are apparent in Fig. 
[2j First, in the low density limit, n e is no longer a linear 
function of /i and, second, the ground state energy [Iq is 
even more lowered with respect to the Eg = case. This 
latter feature indicates that, for fixed A, a single electron 
is more strongly coupled to phonons as Eq increases. 

A more quantitative estimation of the role of SO cou- 
pling on the el-ph properties is given by the electron ef- 
fective mass enhancement to* Jm. This quantity can be 
evaluated from 



= 1 71=0 

TO U)„ 



(21) 



provided sufficiently low temperatures are considered. 
We have checked that for T = 0.02cjo the effective mass 
ratio extracted from Eq. (|2ip is in very good accord 
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with the mass enhancement obtained from the real fre- 
quency self-energy (see next section). In Fig. [3] we report 
our results for to* /to as a function of the electron num- 
ber density n e for the same parameter values of Fig. [5] 
For Eq = we obtain the typical trend for a 2D elec- 
tron gas in the non-crossing approximation: m*/m is al- 
most a constant and approximately equal to the Migdal- 
Eliashbcrg result 1 + A for relatively large densities while, 
for n e — > 0, to* jm decreases towards the one electron 
result^ For E = 5u>q, the mass enhancement follows 
the E = case for densities larger than n e ~ 0.2, cor- 
responding to the range of densities for which n e is pro- 
portional to \i (see Fig. [2]). Instead, for lower values 
of n e , to* /to increases up to a maximum and eventually 
decreases again as n e ^ 0. Higher values of Eq empha- 
size the same trend, with higher and broader maxima of 
to* /to as Eq increases. 

The results plotted in Fig. [3] clearly show how the un- 
derlying diverging DOS, Eq. (j2TJ|) , for E ^ is responsi- 
ble for the enhancement of the effective mass. By reading 
off from Fig. [5] the values of fi corresponding to the den- 
sity values for which to* /to deviates from 1 + A, it is easy 
to realize that the enhancement of to* /to starts when fi 
becomes lower than ~ Eo, that is when the (bare) DOS 
diverges as ^/Eq/E. In this situation, the coupling to 
the phonons is no longer parametrized by A alone, but 
rather by an effective coupling which takes into account 
the strongly varying DOS at low energies^ As a matter 
of fact, for small A, by enhancing Eq the system crosses 
over from a weak to a strong coupling regime, where the 
mass enhancement can be considerably larger than unity. 
It becomes therefore natural to consider signatures of 
such SO induced strong el-ph coupling regime also in 
the spectral properties of the electrons, which can pro- 
vide valuable information testable by tunneling and/or 
photoemission experiments i 23 i 24 ' 26 

IV. SPECTRAL PROPERTIES 

The self-energy for real frequencies could be obtained 
directly from analytical continuation on the real axis of 
Eqs. Ull); (HSJ) and (TTo| . However, since convergence 
is faster on the imaginary axis, in this paper we opt for 
the more efficient method of analytical continuation for- 
mulated in Ref . [3l| . Hence, once T,(iu) n ) has been deter- 
mined from the imaginary axis equations (|14p ~ (|16p . the 
retarded self-energy T,r(uj) = T,(ui + iS) follows from 

S fl (w) = -TX^Djuj - iui m )g(iui m ) 

m 

+^YK w o) + /(^0 - - ^o) 

I n ( w o) + /(^o + u))]g R (u + ui ), 

(22) 

where gn{uj ± ujq) = g{u> ± lu + iS) and n(x) and f(x) 
are the distribution functions for bosons and fermions, 
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FIG. 4: (color online). Real and imaginary parts of the elec- 
tron self-energy for A = 0.5 and electron density n e = 0.1. 
The SO energy is Eq — in (a) and Eq = 5wo in (b). 



respectively. The real and imaginary parts of Sfl(w) are 
plotted in Fig. [4] for A = 0.5, n e = 0.1, E = [Fig. 
@Ja)], and E = 5ojq [Fig. BJb)]. The mass enhancement 
extracted from m*/m = lim w ^o[l — dReT, R (u}) / dco] is 
1.45 for Eo = and 1.72 for Eq = 5ujq, which agree 
with the to* /to values plotted in Fig. [3l For Eq = 0, 
the self-energy displays features typical of the Holstein 
model for a 2D system in the non-crossing approxima- 
tion. Namely, ImS]/?(a;) = for \uj\ < luq, while at larger 
frequencies IuiT,r(lu) ~ — 7rAwo/2. The rapid decrease 
of |ImSfl(w)| at negative frequencies stems from the bot- 
tom band edge. For Eq = 5luq [Fig. [D(b)] the structure of 
£fl(<x>) is more intricate due to the strong energy depen- 
dence of the underlying bare DOS. In fact, for n e = 0.1, 
the value of the (bare) chemical potential fi is well be- 
low Eq = 5uj (see Fig. [5]), and the w-dependence of 
T,r(uj) becomes strongly influenced by the square-root 
divergence of the DOS. This is particularly clear in Fig. 
HJb), where ImE(a;) reproduces for u> < the low-energy 
profile of the DOS shifted by multiples of ojq. This fea- 
ture is characteristic of a strongly-coupled el-ph system 
and is fully consistent with the high value of the mass 
enhancement (to*/to ~ 1.72) for this particular case. 

A global view of the behavior for several values of the 
electron number density and of the SO energy is given in 
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FIG. 5: (color online). Reduced DOS for A = 0.5 (solid lines) and A = (dashed lines). The vertical dotted line at u — 
indicates the Fermi level. The cutoff energy is E c = lOOt^o and the temperature is T = 0.02ajo for the interacting cases (T = 
for A = 0). (a): E = 0, (b): E = 5w , (c): E = 10w 0) (d): E = 20u . 



Fig. HJwhere the reduced DOS for the interacting system 
p(w) = Img R (uj) (23) 

is plotted for fixed A = 0.5. For comparison, we report 
also the bare DOS po(u>), Eq. (|20|) . for the corresponding 
values of n e and Eq. For Eq = 0, Fig. EJa), reduc- 
ing the electron density merely shifts the Fermi level for 
the interacting electron (vertical dotted line) towards the 
bottom of the band. For \u>\ < ujq p(u>) coincides with 
the bare reduced DOS po(uj) = 1 because, as also shown 
in Fig. 21 the imaginary part of the self-energy is zero 
in that frequency range. Compared to the A = case, 
whose DOS has a finite step at the bottom of the band, 
the profile of p(u>) is smeared by the el-ph interaction. A 
similar feature is obtained also for Eq = 5u>q [Fig. E{b)] 
and n e = 0.3 where, now, the square-root divergence of 
po(u>) is rounded-off in p{u>) due to the finite lifetime for 
A = 0.5. However, contrary to the Eq = case, reducing 
n e does not translate to a (more or less) rigid shift of 
the Fermi level but, rather, creates new structures whose 
intensity increases as the Fermi level moves deeper into 
the square-root singularity of Po(lj). This is even more 
pronounced for Eq = IOloq and Eq = 20cjo plotted re- 
spectively in Figs. [5] (c) and (d). For the latter cases, 
the profile of p(u>) for n e = 0.1 is characterized by well 
defined peaks separated by multiples of loq, and whose 
widths decrease as Eq/loq is enhanced. 



Such strong-coupling features are in principle directly 
observable by means of low-temperature tunneling or 
photocmission measurements provided, however, that 
other interactions do not alter significantly the profile 
of p(u>). This could be not the case for example when 
disorder effects arc taken into account, since these tend 
to smear all sharp features of the DOS even at zero tem- 
perature. To investigate this point, we have considered a 
coupling to a short-range impurity potential of the form 
V(r) = Vimp J2i ^( r — R-i)) where are the random po- 
sitions of impurity scatterers. Within the self-consistent 
Born approximation, the resulting self-energy is hence 
given by Eq. (|15p with the impurity term rj mp g(iu; n )/7r 
added in the right-hand side, and Eq. (|22"]) modified ac- 
cordingly. The parameter r imp = l/2r imp = irniV^No 
is the usual scattering rate for zero SO coupling and for 
density of impurities. In Fig. [S]we report the calcu- 
lated reduced DOS p(u>) for A = 0.5, Eq = 10u> , n e = 0.1 
and for several values of r; mp . Also plotted by dashed 
lines for r; mp ^ are the corresponding DOS curves in 
the absence of el-ph interaction (A = 0). Compared to 
the clean limit Fj mp = (top panel), for rather weak 
disorder (r imp = 0.02-E ) the phonon peaks at lu < 
are considerably less sharp and slighly shifted at more 
negative frequencies, but nevertheless still clearly dis- 
cernible. Larger values of Fj mp increasingly smear the 
phonon structures, and gradually the peaks disappear. 
For r imp = 0.2i?o (which corresponds to r imp = 2uj ) the 
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FIG. 6: (color online). Reduced DOS for A = 0.5 (solid lines) 
and A = (dashed lines) for several values of the impurity 
scattering rate Fimp. The temperature is T = 0.02tJo and 
Eq — IOluq, Tie = 0.1 for all cases. 



resulting DOS is basically dominated by the impurity in- 
teraction, and does not deviate much from the A = 
case. 



V. DISCUSSION 

In this section we review the meaning of the approxi- 
mations used in the present work and discuss alternative 
models for the study of the el-ph interaction in strong 
Rashba SO systems. Let us start by considering the lim- 
itations of the non-crossing approximation for the elec- 
tron self-energy. For zero SO interaction, or for Fermi 
energies sufficiently larger than Eq, this is a rather good 
description of the el-ph problem provided A is sufficiently 
small. However, as we have seen, for non-zero SO interac- 
tion the electrons behave as effectively strongly coupled 
to the phonons when the Fermi level is below Eq . This is 
because of the square-root singularity in the DOS when 
Eq 7^ 0. In the Ep < Eq regime, therefore, the non- 
crossing approximation, although making evident the 



trend towards strong coupling, may not be adequate for 
a more trustworthy description of the system. It is in- 
structive at this point to consider the limiting situation 
where only one electron is present. By using second-order 
perturbation theory (that is simply the non-crossing ap- 
proximation with the bare electron Green's function) it 
is easy to evaluate the mass enhancement factor. At zero 
temperature, and setting for simplicity E c /Eq = oo, this 
is given by 



m* A A [Eq~ ( [Eq' 
= 1 + — + — \ / — arctan \ / — 



(24) 



It is clear from the above expression that for Eq 3> ojq the 
mass enhancement factor is governed by an effective cou- 
pling, say A, proportional to Xy 1 Eq/ujq, amplified with 
respect to A by the square-root divergent DOS. Equation 
(|2"l)l clarifies also that the relevant adiabatic parameter 
for Ep <C Eq is ujq/Eq, rather than uq/E c (where E c 
plays the role of the bandwidth), and that the effective 
coupling A increases as luq/Eq — > 0. Consequently, in the 
adiabatic limit luq/Eq = perturbation theory breaks 
down for any finite A because A = oo. This leads us 
to suspect that, in analogy with the adiabatic limit of 
the one-dimensional lattice Holstein model^ the ground 
state of a single electron for luq/Eq = is always a bound 
polaron. However, for ujq/Eq ^ 0, the effective coupling 
A ~ Eq/ojq is finite, which permits to estimate a 
rough range of validity of Eq. ([M|) . Indeed, contributions 



of higher orders of perturbation theory become negligi- 
ble as long as A <C 1, corresponding to A <C \Jloq/Eq, 
consistent with the results on the one-dimensional Hol- 
stein model of Refs. 33li34j . which show better agreement 
between perturbation theory and exact diagonalization 
results as loq increases. 

Also for low but finite electron densities it is possi- 
ble to interpret the coupling to the phonons in terms 
of an effective el-ph coupling A which grows as Eq in- 
creases. For example, in the range of electron densities 
n e = 0.1-0.2, the mass enhancement factor plotted in 
Fig. for A = 0.5 may be interpreted by an effective 
Migdal-Eliashbcrg formula m* jm = 1 + A, where Awl 
for Eq = lQujQ and A w 1.3 for Eq = 20wo- However, 
contrary to the one-electron case discussed above, now 
A depends on the Fermi energy Ep. Indeed, provided 
that wo < Ep < Eq, the effective coupling turns out 
to be of order A ~ \\J Eq/ Ep, where the square-root 
term stems from the singularity of the DOS, Eq. (f20"|) . in 
analogy with the general definition of the effective el-ph 
interaction in the presence of a van Hove singularity. 35 
At this point it is possible to estimate the validity of the 
self-consistent non-crossing approximation for the self- 
energy considered in the previous sections. In fact, ac- 
cording to Migdal's theorem generalized to systems with 
diverging DOS,— the el-ph vertex correction factors be- 
yond the non-crossing approximation are at least of order 
P = Xujq/Ep, so that neglecting them would introduce 
an error of order P. Estimates of P for the different 
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cases discussed in this paper can be obtained by evalu- 
ating A ~ 1 — m*/m from FigJ3] In this way the Fermi 
energy is roughly given by Ep ~ (\/\) 2 Eq, which then 
can be inserted in the definition of P. For the low density 
value n e = 0.1 we find P > 1 for Eq = 20ujq and P ~ 0.5 
for Eq = 10wo, showing that the non-crossing approxima- 
tion is quantitatively inaccurate in this case. However, 
already for n e = 0.2, for which effectively strong-coupling 
features are apparent from FigJ3] and FigdJ the contri- 
butions of the vertex corrections drop to P ~ 0.4 and 
P ~ 0.2 for Eq = 20wo and Eq = lOwoi respectively. In 
this situation, the non-crossing approximation is fairly 
reliable and its accuracy improves as n e is further en- 
hanced and/or Eq/loq is reduced. 

Let us turn now to discuss the general form of the self- 
energy for the case in which the el-ph matrix element is 
momentum dependent. Here we consider the situation 
in which the momentum dependence is only through the 
modulus q of the momentum transfer q, as is the case, 
for example, with the 2D Frohlich model, for which the 
coupling goes like 1/ ' -Jq~. As shown in the Appendix, 
a fully general expression of the self-energy valid also 
beyond the non-crossing approximation is: 

S(k, iw n ) = T,i(k,iu) n )l + S 2 (fc, iw„)f2 k ■ cr, (25) 



where Ei and £ 2 are scalars. Compared to Eq. (|T5|) . the 
above expression has an additional term which is off- 
diagonal in the spin subspace, renormalizing therefore 
the SO coupling. This term disappears (£ 2 = 0) only 
when the el-ph matrix element is momentum indepen- 
dent, like in the Holstein model, and at the same time 
the self-energy is evaluated in the non-crossing approxi- 
mation. In all other cases, like, e.g., the Frohlich model 
in the non-crossing approximation, E 2 is nonzero. For 
sufficiently large values of Ep, such that the weak SO 
limit Eo/Ep <C 1 holds true so that the Fermi level lies 
far above the ID singularity of t he DOS , £ 2 turns out to 
be of order AwA so / Ep cx Aa>o-\/ Eq/Ep, and can be dis- 
regarded in comparison with Si « Au>o- On the contrary, 
when Ep/Eo < 1, Si and S 2 have comparable magni- 
tude, and the full momentum and frequency dependent 
of both terms must be considered for a consistent evalu- 
ation of the el-ph effects. 



VI. CONCLUSIONS 

In this paper we have addressed the role of the Rashba 
SO interaction in the properties of a coupled el-ph gas in 
two dimensions. By using a self-consistent non-crossing 
approximation for the electron self-energy, we have stud- 
ied the mass enhancement factor and the spectral prop- 
erties. We have shown that, for sufficiently strong SO 
interaction, the electron becomes strongly coupled to the 
phonons even if el-ph coupling A can be classified as weak. 
We identify this behavior as being due to a topological 
change of the Fermi surface for strong SO interaction, 
which gives rise to a square-root singularity in the DOS 



at low energies. Signatures of such effectively strong el- 
ph coupling are found in the mass enhancement factor, 
which becomes as large as m*/m ~ 2 for el-ph coupling 
of only A = 0.5, and in the energy dependence of the 
interacting DOS, displaying low energy peaks separated 
by multiples of the phonon energy loq. This latter feature 
could be tested experimentally by tunneling or photoc- 
mission experiments in systems where the Fermi level 
can be tuned to approach the square-root singularity of 
the DOS. We have then discussed limitations of the non- 
crossing approximation approach and possible generaliza- 
tions of the theory for momentum-dependent el-ph ma- 
trix elements. Since the problem of el-ph coupling in the 
presence of SO interaction is relevant for several systems 
such as metal and semimetal surface states, surface su- 
perconductors, or low-dimensional hcterostructures, and 
given the current interest in spintronic physics, we hope 
that our work will stimulate further investigations. 



APPENDIX 

In this appendix we evaluate the form of the elec- 
tron self-energy when the el-ph interaction is momentum 
dependent. In particular, we consider the hamiltonian 
H = Ha + H p h, where Hq is the Rashba spin-orbit hamil- 
tonian of Eq. JT]) and 



H 



ph 



E 



uj a, 'a 



E ffk-k'cj CQ c k - Q (aj c _ k , +a k /_ k ), 



kk', 



(A.l) 

where g q is the el-ph matrix element which we assume 
depends only on the modulus of momentum transfer q. It 
is convenient to rewrite H in terms of the eigenvectors of 
Hq, whose annihilation operators d ks (s = ±) are related 
to c ka through 



Cki 



1 



1 -ie-^f 



V2 V 1 ie 



C k | 



(A.2) 

where tp is the azimuthal angle of k. In this basis, Hq is 
diagonal with dispersion relation given by Eq.([3]) while 
Hph becomes 

Hph = y^u>oalag+ ^ M k-k' d w s d ^'^'( a k-k' + fl k'-k), 

(A.3) 



kk'ss' 



with 



2 k _ k , = 5k-k' 



1 + ss'e^-*'") 



(A.4) 



By applying Wick's theorem, it turns out that, to all 
orders of the el-ph interaction, the Green's function in 
the chiral basis has zero off-diagonal components so that, 
if G±(k,r) = -(T T d k ±(r)4±(0)), the matrix Green's 
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function in the original spin sub-space becomes 



G(k, iuj n ) = T 



t /G+(k,icj n ) 



GL(k,icj n ) 1 
= \ V(l + «n k -o-)G.(k,iw n ). (A.5) 

Consequently, by using Dyson's equation ([9]), the matrix 
self-energy in the spin sub-space is 

S(k, iu n ) = Si(k, iu n )l + £ 2 (k, iw n )fi k • <x, (A.6) 

where 

v /, ■ s S+(k, zo; n ) + (-)£_(k, w„) 

S 1(2) (k,«o; n ) = . (A.7) 



In order to obtain Eq. (|25[) . it suffices to demonstrate that 
the momentum dependence of self-energy in the chiral ba- 
sis, £±(k, iio n ), is only via k — |k|. This is accomplished 
by noticing that the el-ph matrix element in the chiral 
basis, Eq. (|A.4[) . depends on the direction of the momen- 
tum transfer k — k' solely through <p — ip'. Hence, if the 
electronic dispersion depends only on the modulus of the 
momentum, as is the case with Eq. , a general self- 
energy diagram in the chiral basis will be independent 
of the direction of k which, by using (|A.7|) . is consistent 
with Eq. ([25]) . 
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